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ABSTBRCr 


This  report  deals  with  design  optimization  of  Internal- 
Expanding  Rio  Brakes.  A  computer  program  was  developed  to 
calculate  the  actuating  force,  torque,  stopping  time  and 
drum  temperature.  The  drum  temperature  is  calculated  by  the 
finite  difference  method. 

A  comparison  of  results  has  been  made  using  a  simplified 
equation  that  is  in  common  use  in  engineering  texts. 


Numerical  optimization  is  shown  to  be  a  convenient  tool 
for  brake  design. 
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A.  ENGLISH  LETTER  SYMBOLS 

a  Distance  from  pivot  to  the  center  of  rotation  (m) . 

A  Area  of  one  lining  shoe  (m^). 
b  Width  of  friction  material  (m)  . 

B^  Biot  modulus, 
c  Specific  heat  (J/Kg“°C). 

C  Thermal  capacity  (J/®C) . 

d  Distance  from  actuating  force  to  the  hinged  pin  (m)  . 
dc  Rata  of  deceleration  (m/sec!  )  . 

E  Kinetic  energy  (J)  . 
f  Frictional  force  (N) 

P  Actuating  force  (N)  . 

Pfl  Fourier  modulus, 
g  Gravity  constant  (m/seci  ) . 

h  Convection  heat  transfer  coefficient  (W/a^-C). 

k,  Thermal  conductivity  (i/a-'’C)  . 

Friction  moment  (N-m)  . 

M  Normal  moment  (N-m)  . 

n 

N  Normal  force  (N)  . 

p  Pressure  between  lining  and  drum  at  any  point  (N/m^). 
p  naximum  pressure  between  lining  and  drum  (N/m^). 

Q  Heat  generated  (W) . 

r  Inside  drum  radius  (m) . 

R  Wheel  radius  (m)  . 

R^^  Thermal  resistance  (°C/W). 

t  Time  (sec.) 

tk  Thickness  (m) . 

T  Temperature  (°C)  . 

To  Torgue  (N-m)  . 

V  Velocity  (m/sec.). 

Vq  Volume  (mM  • 

W  Vehicle  weight  (N)  . 
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B.  uoTirroN 

R.  .  The  thermal  resistance  between  node  i  and  the 
ij 

adjoining  node  j. 

The  temperature  of  node  i  at  time  step  p. 


C.  GREEK  LETTER  SYMBOLS 

9  The  angle  between  the  hinged  pin  and  an  element 
area  on  the  lining. 

9  The  angle  at  which  the  pressure  between  the  lining 

3, 

and  drum  is  maximum. 

U  Friction  coefficient. 

u  Cold  friction  coefficient, 
c 

Hot  friction  coefficient, 
a  Thermal  diffusivity  (m^/sec.). 

A  Finite  increment. 
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1.  IHTaODnCTION 

Brakes  are  aechanical  devices  for  retarding  the  motion 
of  a  vehicle  or  machine  by  means  of  friction.  Because  of  the 
similarity  of  their  functions,  many  clutches  may  also  be 
included  here,  assuming  centrifugal  forces  are  accounted 
for. 

A  simplified  dynamic  representation  of  a  brake  is  shown 
in  Pig.  1.  Two  masses  with  inertias,  Ij  and  ,  rotating 
at  the  respective  angular  velocities  oji  and  6U2  (one  of  which 
may  be  zero) ,  are  to  be  brought  to  the  same  speed  by 
engaging  the  brake. 

The  friction  brake  has  three  basic  elements;  two 
opposing  friction  surfaces  and  a  mechanism  for  forcing  the 
friction  surfaces  into  contact.  Whenever  a  friction  brake  is 
engaged  to  join  two  members  having  relative  motion,  there  is 
a  period  of  slip  which  may  last  several  seconds.  This  slip 
is  one  of  the  chief  merits  of  the  friction  brake;  it  absorbs 
saocks  and  prevents  excessive  torsional  stresses  on  the 
power  transmission  system.  On  the  other  hand,  slip  is  the 
limiting  factor  in  friction  clutch  and  brake  performance; 
for  heat  is  generated  in  proportion  to  slip,  torque 
transmitted,  and  period  of  slip. 

The  following  parameters  are  of  interest  in  analyzing 
the  performance  of  these  devices; 

1.  The  actuating  force. 

2.  The  torque  transmitted. 

3.  The  temperature  rise. 

4.  The  slip  time. 

This  report  deals  with  Internal-Expanding  Rim  Brakes, 
This  formulation  also  applies  to  internal-expanding  clutches 
if  centrifugal  forces  are  accounted  for. 
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II.  INTSRSAL-BXPANDING  3I?1  CLaTCHES  AND  BRAKES 

A.  GENERAL  HACHANICAL  PRINCIPALS 

A  bralce  or  clutch  assembly,  uses  a  brake  shoe  to  which 
is  attached  a  frictioa  material,  called  lining.  The  lining 
is  riveted  or  bonded  to  the  brake  shoe  as  shown  in  Fig.  2. 
The  brake  shoe  is  pivoted  at  a  fixed  point  and  the  other  end 
is  subjected  to  a  force  which  presses  the  shoe  in  contact 
with  the  drum.  The  force  between  the  brake  and  the  drum  is 
radial  as  the  drum  rotates.  If  a  point  on  the  rotating  drum 
surface  first  makes  contact  with  the  shoe  at  the  end  nearest 
the  pivot,  the  shoe  is  termed  a  "trailing  shoe".  If  it  first 
makes  contact  at  the  other  end  the  shoe  is  tamed  "leading 
shoe",  tha  latter  giving  a  higher  braking  torque  than  the 
former  for  a  given  braking  force. 

The  friction  batween  the  lining  and  the  drum  creates 
haat  which  is  basically  the  conversion  of  energy  of  motion 
of  the  vehicle  or  machine  to  thermal  energy  at  the  friction 
surfaces,  namely  the  lining  and  the  drum.  This  heat  is  then 
dissipated  and  absorbed  by  the  drum  by  conduction, 
convection  and  radiation  into  the  atmosphere. 

B.  FRICTION  PaNDAMESTALS  AND  NATEEIALS 

Friction  mechanisms,  such  as  brakes,  are  systems  for 
converting  mechanical  energy  into  heat.  Several  basic 
factors  affect  friction  and  wear  of  materials  used  in  brake 
systems.  The  main  factors  are  temperature,  pressure,  speed, 
surface  roughness,  and  type  of  material.  Some  organic  or 
molded  friction  materials  show  no  change  in  friction 
characteristics  with  pressure,  while  others  such  as 
sinter ed-metal  materials  decrease  in  friction  coefficient  as 
pressure  is  increased.  For  metallic  friction  materials  there 
is  also  a  decrease  in  coefficient  of  friction  as  speed 
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increases.  Temperatare  effects  upon  the  coefficient  of 
friction  vary  widely  with  the  type  of  materials  used. 

In  a  two-shoe  internal  expanding  brake  there  is  a 
tendency  for  the  brake  drum  to  deform  under  hard 
application.  Drums  become  elliptical  and  the  force  to  do 
this  is  guite  high  and  contributes  to  friction  force. 

A  brake  or  clutch  friction  material  should  have  the 
following  characteristics  to  a  degree  which  is  dependent 
upon  the  severity  of  the  service: 

1.  A  high  and  uniform  coefficient  of  friction. 

2.  The  ability  to  withstand  high  temperarures,  together 
with  good  heat  conductivity. 

3.  Properties  which  are  not  affected  by  environmental 
conditions  such  as  moisture. 

U.  Good  resiliency. 

5.  High  resistance  to  wear,  scoring  and  galling, 

C.  BRAKE  ORONS 

One  of  the  primary  functions  of  a  brake  drum  is  that  of 
absorbing  and  dissipating  the  heat  developed  during  the 
application  of  the  brake.  A  brake  drum  is  a  heat  sink  into 
which  heat  goes  after  it  is  created  by  the  rubbing  friction 
of  the  brake  lining  contact  to  drum.  The  brake  shoe  and 
lining  permanently  fixed  on  the  axle,  when  actuated, 
contacts  the  drum  under  pressure  to  cause  the  friction  to 
stop  the  vehicle.  The  energy  of  motion  of  a  vehicle  is 
converted  to  thermal  energy  b)  the  brake  assemblies.  A  brake 
drum  must  have  the  capacity  to  absorb  and  dissipate  this 
heat  energy  within  the  limits  of  the  brake  heat  input.  If 
this  is  not  the  case,  the  drum  expands  and  the  brakes  fade 
or  fail.  The  greater  the  mass  of  the  drum,  the  more  heat  it 
can  absorb  and  store  until  such  time  as  the  heat  can  be 
dissipated  by  convection  and  radiation  [Ref.  1]. 

An  ideal  brake  drum  would  have  the  following 

characteristics; 
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1.  Higli  structural  strength  to  resist  bursting  forces. 

2.  Oniform  coefficient  of  friction. 

3.  Hard  surface  to  resist  scoring. 

4.  High  heat  conductivity  to  rapidly  conduct  heat  away 
from  braking  surfaces. 

5.  High  emissivity  factor  to  radiate  heat  from  the  drum 
surface  to  the  atmosphere. 

6.  High  heat  storage  capacity  to  store  heat  from 
successive  brake  applications  until  it  can  be 
dissipated. 

7.  Good  aachinability  to  permit  boring  of  the  drum. 

D.  STATIC  AND  DYNAMIC  ANALYSIS 

1 .  Assumptions 

In  developing  the  equations,  the  following 
assumptions  have  been  made; 

a.  The  pressure  at  any  point  on  the  shoe  is  proportional 
to  the  moment  arm  of  this  point  from  the  pivot. 

b.  The  affect  of  centrifugal  force  may  be  neglected. 

c.  The  shoe  is  assumed  to  be  rigid. 

d.  The  friction  coefficient  is  a  linear  function  of 
temperature  and  it  does  not  vary  with  pressure,  wear 
and  environment. 

2.  Pressure  Concept 

To  analyze  an  internal  shoe  refer  to  Fig.  2,  which 
shows  a  shoe  pivotal  at  a  fixed  point  with  the  actuating 
force  acting  at  the  other  end  of  the  shoe.  The  mechanical 
arrangement  does  not  permit  pressure  to  be  applied  at  the 
pivot,  therefore  the  pressure  at  this  point  is  zero.  If  the 
shoe  rotates  through  a  small  angle  about  A,  the  radial 
movement  of  any  point  on  the  arc  of  contact,  is  proportional 
to  the  moment  arm  of  this  point  from  the  pivot.  Assuming 
that  the  material  of  the  brake  lining  and  support  obey 
Hooke's  law,  the  pressure  at  this  point  will  also  be 
proportional  to  this  moment  arm.  The  distance  is 
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proportional  to  sin  9  .  Therefora,  the  relations  between 
pressure  at  any  point  and  the  maxiniuin  pressure,  p  ,  will  be 
given  by  the  following  formula; 

P  P 

sin9  sin9 

a 

Prom  this  formula  it  can  be  seen  that  the  frictional 
material  at  the  heel,  contributes  very  little  to  the  braking 
action,  therefore  it  is  better  to  begin  the  friction 
material  at  an  angle  9i  greater  than,  say  0.15  rai.  It  can 
be  seen  also  that  the  pressure  will  be  marimum  when  9  =90® 
or  if  the  toe  angle  92  is  less  than  90  then  the  pressure 
will  be  maximum  at  the  toe.  For  good  performance  it  is 
recommended  to  concentrate  as  much  frictional  material  as 
possible  in  the  neighborhood  of  the  point  of  maximum 
pressure  [Ref.  2]. 

3.  kctuating  Force  and  Torgue  Calculation 

Prom  Fig.  2,  it  can  be  seen  that  the  differential 
normal  force  on  an  element  area  of  the  lining  will  be; 

dN  =  pdA  (2) 

where  dA  is  an  area  element  of  the  lining  and  it's  magnitude 
is; 


dA  =  rbd9 


(3) 


In  Equation  3,  r  is  the  inside  drum  radius  and  b  is  the  drum 
width.  Substituting  for  p  and  dA  gives; 


p  brsin9 

dN  =  ■;  - —  d9 

sin9 


(4) 


At  the  same  point  the  differential  frictional  force  is; 


df  =  udN 

where  u  is  the  coefficient  of  friction. 


(5) 
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The  actuating  force,  F,  can  be  calculated  using  the  fact 
that  the  suanation  of  the  aoments  about  the  hinge  pin  is 
zero.  The  moment  due  to  frictional  forces  is; 

M-  =fl;®2(r-cos9)df  (6) 

I  w  1 


where  a 
rotation 
from  0j 


is  the  distance  from  the  pivot  to  the 
.  Substituting  the  value  of  dF  and  i 
to  02  gives; 

a 

Q — { (cos0  I~cos0  2)+-5p(sin^0  1-Sin^0  2  )  } 


center  of 
ntegrating 


(7) 


a 

where  y  is  assumed  to  be  constant  along  the  lining. 
Similarly  the  moment  due  to  normal  forces  is  given  by; 


=„/®^asin0dN 
n  □  1 


(8) 


Substituting  the  value  of  dN  and  integrating 
from  01  to  02  gives; 
p  bra 

M  =  V— A-{0. 5(0 2-01 )-O.25(sin9  2-sin0i)  (9) 

n  sino_ 
a 

rhe  actuating  force  must  balance  the  moments, therefore; 


F  a: 


M  -M„ 
n  f 


(10) 


where  d  is  the  distance  from  the  hinge  to  the  point  of 
application  of  F.  The  torgue  applied  to  the  drum  by  the 
brake  shoe  is; 


To  =9{®"rdf  (11) 

After  substituting  the  value  of  df  and  integrating  ; 

UP„br2 

To  =  — (cos9i-cose2)  (12) 
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4.  Rate  of  Heat  Senerated  ani  Deceleration  Calculation 
The  differential  rate  of  heat  generated  by  an 
eleoent  area  of  the  lining  is  equal  to  the  velocity  of  the 
inside  surface  of  the  drua  relative  to  the  lining  ,  times 
the  differential  frictional  force  acting  on  the  element 
area; 


dQ  =  V^df  (13) 

Assuaing  the  brake  is  on  a  vehicle  wheel  with  a  radius  of  R, 
the  inside  surface  velocity  is  equal  to; 


where  V  is  the  velocity  of  the  vehicle  and  is  a  function  of 
time. 

If  V=V  (t)  then  7  =V  (t)  and  the  heat  generated  will  be  also 

r  r 

a  function  of  time.  Substituting  the  values  of  and  df  and 
integrating  from  9 i  to  02  #  we  get  the  following  formula  for 
the  heat  generated  at  any  time  t, 

Q(t)  =  ■gj^^(|')(cos9i-cose2)V(t)  (15) 

a 

The  kinetic  energy  of  a  vehicle  of  weight  W  is  given  by; 

E  =  |(|)V"  (16) 

dote  that  if  the  brake  is  on  a  four  wheel  vehicle,  there 
will  be  eight  shoes.  Assuming  all  are  leading  shoes,  each 
will  stop  one-eight  of  the  vehicle  weight,  so  W/8  must  be 
used  in  Equation  (1d|  .  The  rate  of  change  in  the  kinetic 
energy  is; 


dt 


From  the  energy  conservation  law  the  rate  of 
kinetic  energy  is  equal  to  the  heat  generated; 


(17) 

change  in  the 
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(18) 


Q(t) 


dt 


Substituting  the  value  of  Q  (t)  aad  dE/dt,  it  is  seen  that 
the  velocity  V(t)  cancels  and  so  the  deceleration  is  not  a 
function  of  tiae.  Therefore  the  deceleration,  dc,  is; 

Hxr  P  by  2 

a 

The  velocity  at  any  riae  is; 


V  =  V^-dct 


(20) 


where  is  the  initial  velocity.  Substituting  the  velocity 
i.i  Equation  (16)  ,  yields  the  rate  of  heat  generated  as  a 
function  of  time. 


P„bu  _2 

Q(t)  *  ^~^(J)(COS0i-COS02)(V^-dCt) 


(21) 


In  this  study  the  friction  coefficient  was  taken  as  constant 
up  to  a  temperature  of  90®C  and  after  90‘’c,  decreases 
linearly  to  zero  at  a  specified  temperature, 

u  T  <  90  °C 


P  =  ) 


‘‘c-  -^<^-90) 


90  °C  <  T  <  T 

—  m&x  • 


(22) 


T  >  T 


max . 


where  is  the  cold  coefficient  of  friction  and  y^j  is  the 
hot  coefficient  of  friction. 


E.  SORF&CE  TEnPERATOaS  CALCOLATION 

Since  the  function  of  a  brake  is  to  convert  kinetic 
energy  into  heat,  surface  temperatures  of  brake  linings  and 
drums  are  most  important.  Therefore  it  is  necessary  to  know 
the  temperature  of  the  mechanism  during  and  after  any  stop. 
The  temperatures  were  calculated  by  the  finite  difference 
method. 
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.  Assumptions 

a.  One  dimensional  heat  flow-Ihe  heat  flow  is  from  the 
inner  surface  to  the  outer  surface  of  the  drum. 

b.  Constant  heat  transfer  coefficient. 

c.  No  heat  dissipated  by  radiation. 

d.  The  heat  is  generated  on  the  inner  surface. 

2.  Temperature  Analysis 

a.  Theory 

The  differential  equation  to  be  solved  in  order 
to  find  the  temperature  in  the  drum,  based  on  the 
assumptions,  is; 


3^T  .  Q  _  ,lx3T 
3x''  k  ~  ^a^3t 


(23) 


with  the  following  boundary  conditions: 
at  x-0  heat  is  generated, 

at  xstk  heat  is  transfered  to  the  atmosphere  by 
convection. 

In  the  equation  above  k  is  the  thermal  conductivity,  a  is 
the  thermal  diffusivity,  t  is  time  and  th  is  the  drum 
thickness.  This  eguation  can  be  solved  by  the  finite 
difference  method  'aef.  3].  The  finite  difference  model 
used  here  is  shown  in  Pig.  3.  The  rate  of  change  with  time 
of  the  internal  energy  of  a  node  i  is  approximated  by; 
,j.P+1_tP 

§  =  OCAV.^^  (24) 


where  p  is  the  density,  c  is  the  specific  heat  and  is 
the  drum  volume. 

Now  define  the  thermal  capacity  as 


Ci  =  p^c^AVo.  (25) 

The  forward  difference  equation  for  all  nodes  and  boundary 
conditions  is; 
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(26) 


»jiP_fpp 


P+1_tP 


th,  ij 


C.- 

1 


'j't' 

i 


i 


At 


where  .  .  is  the  theraial  resistance 

n+1 

Solving  the  above  equation  for  gives; 

th.ij  1  1  th,ij 


(27) 


The  thermal  resistance  can  be  calculated  from  the  geometry 
and  boundary  conditions  'Ref.  3].  To  ensure 

stability  At  must  be  equal  or  less  than  the  following 
nodal  relation; 


At  <  ( 


C. 

1 


y _ L 


■) 


(28) 


th,  ij 

With  the  assumptions  made,  the  drum  can  be  viewed  as  an 
infinite  plate,  with  heat  generated  at  the  surface  of  the 
first  node,  as  shown  in  Pig.  3.  It  is  assumed  that  in  every 
drum,  there  are  two  shoes  and  that  both  are  leading  shoes. 
Therefore,  two  times  must  be  taken. 


nP+1 


fjtP 


i  =  (2QP+I5J - t  (1-  - )tP 

^  *^th,ij  ^i  ^i  ‘^th,ij  ^ 

Formulation 


(29) 


0. 


In  the  computer  program  5  nodes  were  taken.  In 
order  to  check  accuracy,  the  program  was  run  with  7  and  10 
nodes.  In  each  case  the  result  was  the  same  within  5  “c. 
The  heat  is  generated  in  the  inner  drum  surface.  Therefore  Q 
appears  in  the  formula  of  temperature  in  the  first  node  and 
for  all  the  other  nodes  Q  is  equal  zero.  With  the 
assumptions  mentioned  above,  the  heat  transfer  through  the 
drum  is  solved  as  a  heat  transfer  problem  through  an 
infinite  plate,  with  heat  generation  at  the  inner  surface 
and  with  a  heat  convection  boundary  on  the  outer  surface  as 
shown  in  Pig.  3.  Squation  (29)  can  be  simplified  using  two 
dimensionless  parameters,  Biot  and  Fourier  modulii. 
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B.  -  ^ 

(30) 

1  k 

Fo  _  aAt 

(31) 

(Ax)^ 

rhe  final  equations  for  calculating 
nodes  now  become; 

For  the  first  node; 

2Q^At 

—  +  (1-2Fo)tP  +  2FoTP 


For  the  interior  nodes; 


the  temperatures  at  the 


(32) 


For  the  Last  node; 

tP*1  -  2F.{tP_^^3jT.*  -Bi-1)TP)  (34) 


F.  BRAKE  DOTY  CYCLE 

In  addition  to  tae  parameters  mentioned  above  the  design 
of  a  brahe  depends  on  the  initial  speed,  final  speed,  number 
of  stops,  and  the  rest  time  between  each  stop.  In  this 
analysis  a  general  duty  cycle  was  considered  so  that  the 
initial  speed,  final  speed  and  the  acceleration  period 
between  stops  can  be  different  for  each  part  of  the  design. 

In  the  design  examples  presented  here,  a  vehicle  was 
stopped  four  consecutive  times  with  the  following  cycle; 


Initial 

Final 

Rest 

Speed 

Speed 

m/sec. 

m/sec. 

sec. 

1 

25.0 

0,0 

20.0 

2 

25.0 

0,0 

20.0 

3 

25.0 

0.0 

20.0 

4 

25.0 

0.0 

> 

21 


in.  opTiaizMiON 


4.  INTRODUCTIO!! 

Enginsaring  analysis  using  the  digital  computer  has 
become  commonplace.  It  is  less  common  to  use  the  computer 
to  maka  the  actual  design  decisions,  such  as  sizing  of 
structural  members  or  placement  of  mechanical  linkages.  This 
may  be  largely  attributed  to  the  fact  that  fully  automated 
design  requires  techniques  that  are  unfamiliar  to  much  of 
the  engineering  community. 

In  many  enginaering  problems,  it  is  necessary  to 
determine  the  minimum  or  maximum  of  a  function  of  several 
variables,  limited  by  various  linear  and  nonlinear 
inequality  constraints.  It  is  seldom  possible,  in  practical 
applications,  to  solve  these  problems  directly,  and 
iterative  methods  ace  used  to  obtain  the  numerical  solution. 
Machine  calculation  of  this  solution  is,  of  course, 
dasirabla.  The  COHMIN  program  is  available  to  solve  a  wide 
variety  of  such  problems  [Ref.  4]. 

COSMIS  is  a  F0RTR4M  program,  in  subroutine  form,  for  the 
minimization  of  a  multi-variable  function  subject  to  a  set 
of  inequality  constraints.  The  basic  optimization  algorithm 
is  the  Method  of  Feasible  Directions  [Ref.  5].  The  user  must 
provide  a  main  calling  program  and  an  external  routine  to 
evaluate  the  objective  and  constraint  functions  and  to 
provide  gradient  information.  If  analytic  gradients  of  the 
objective  or  constraint  functions  are  not  available,  this 
information  is  calculated  by  finite  difference.  While  the 
program  is  intended  primarily  for  efficient  solution  of 
constrained  problems,  unconstrained  function  minimization 
problems  lay  also  be  solved,  and  the  Conjugate  Direction 
Method  of  Fletcher  and  Reeves  is  used  for  this  purpose 
[Ref.  6]. 
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B.  DEFIHITION  3F  TEHMS 

Host  disciplines  have  a  'iniqua  set  of  nonenclature  ased 
to  describe  the  concepts  within  that  discipline.  Soae  of  the 
coanonly  used  terms  in  numerical  optimization  are  summarized 
here. 

Objective-  The  value  of  the  function  which  is  to  be 
minimized  or  maximized  during  the  optimization  process, 
synonyms  are  cost,  merit  and  payoff.  The  common  mathematical 
designation  is  P(X).  In  the  present  study  the  objective  was 
to  minimize  the  material  in  the  brake  drum. 

Design  variables-  The  parameters  to  be  changed 
during  the  optimization  process  in  order  to  minimize  or 
maximize  the  value  of  the  objective  function.  Synonym; 
decision  variables.  The  common  mathematical  designation  is 
the  vector  X.  Design  variables  considered  in  this  study 
include,  drum  thickness,  width,  the  angle  between  the  hinged 
pin  and  the  end  of  the  lining,  and  the  distance  from  the 
pivot  to  the  center  of  rotation. 

Inequality  constraints-  Dne-sided  conditions  which 
must  be  mathematically  satisfied  for  the  design  to  be 
acceptable.  The  common  mathematical  term  is  G(X)<0  or 
G(X)>0.  If  the  ineqality  condition  is  satisfied  on  G  (X)  ,  the 
design  is  acceptable,  (feasible).  If  it  is  not  satisfied, 
the  design  is  not  acceptable  (infeasible) .  Constraints 
considered  here  include,  vehicle  stopping  time,  maximum  drum 
temperature,  and  actuating  force. 

Side  constraints-  Gpper  and  lower  bounds  on  the 
individual  design  variables  X  .  The  common  mathematical 
representation  is  x}<  X^<  X^. 

Design  space-  The  n-dimensional  mathematical  space 
spanned  by  the  vector  of  design  variables  X. 

Xctive  constraint-  Constraint  G.  (X)  is  called  active 

J 

if  its  value  is  zero  (or  near  zero  for  computational 
purposes) . 
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Inactive  constraint-  Constraint  G.(X)  is  inactive  if 

J 

G  (X)<0. 

J  . 

Tiolated  constraint-  Constraint  G.(X)  is  violated  if 

J 

G  (X)>0. 

j 

C.  THE  OPTIMIZATION  PROCESS 

The  general  design  optimization  problem  can  be  stated 
mathematically  as  follows:  Find  the  set  of  variables  X^, 
i=1,2....n,  which  will 

Minimize  P(X  )  (35) 

Subject  to: 

G.(X)  ■<  0  j  =  l,2 . m  (36) 

•D 

<  X.  <  xY  i=l,2 . n  (37) 

1—1-1 


Vector  X  contains  the  set  of  independent  design 
variables  X^,  i=1,2....n.  X  may  represent,  for  example 
width,  thickness,  and  angles  in  the  brake  optimization.  The 
objective  function  used  here  is  the  drum  volume. 

Equation  (36)  defines  the  ineqality  constraints  imposed 
on  the  design.  For  example,  if  the  temperature  on  the  inner 
drum  surface  must  not  exceed  a  specified  value  T,  the 
associated  design  constraint  becomes,  in  normalized  form 


1  <  0 


(38) 


The  lower  and  upper  bounds  on  the  design  variables, 
given  by  Eq. (37)  ,  limit  the  region  over  which  the  functions 
P(X),  and  G  (X)  are  defined.  These  constraints  are  often 


referred  to  as  side  constraints  because  they  form  the  sides 
or  bounds  of  the  n-dimensional  spaced  spanned  by  the  design 
variables  X. 


If  all  the  inequalities  of  Egns. (36)  and  (37)  are 
satisfied,  the  design  is  said  to  be  feasible;  if  any  of 
these  conditions  are  not  satisfied,  the  design  is  not 
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feasible.  If  F(X)  is  a  ainiatia  and  the  design  is  feasible, 
it  is  also  optimum,  or  at  least,  a  relative  optimum.  Note 
that  because  the  objective  and  constraints  may  be  nonlinear, 
there  may  be  multiple  minima  in  the  design  space  that  cannot 
be  identified  using  current  methods.  While  this  is  a  matter 
for  concern,  since  it  is  desired  to  find  the  true  optimum, 
it  must  be  remembered  that  the  same  mathematical  conditions 
exist  if  the  design  process  is  not  automated.  However,  using 
optimization  techniques,  it  is  a  simple  matter  to  restart 
the  optimization  from  several  initial  points  in  the  design 
space  and  thereby  improve  the  probability  of  obtaining  the 
tcue  optimum  design,  a  process  that  would  be  quite  time- 
consuming  in  manual  design. 

Equations  (35) -(37)  define  the  nonlinear  constrained 
optimization  problem.  If  Egs. (36)  and  (37)  are  not  imposed 
on  the  design,  the  optimization  problem  is  defined  by 
Eg.  (35)  alone  and  is  therefore  an  unconstrained  minimization 
problem. 

Host  aonlinear  optimization  algorithms  update  the  vector 
of  design  variables  by  the  iterative  relationship; 

+  aS*^  (39) 

where  q  is  the  iteration  number,  vector  S  is  the  direction 
of  search  in  the  design  space,  and  the  scalar  a  is  referred 
to  as  a  move  parameter  which, together  with  S,  determines  how 
much  the  vector  X  is  changed  luring  the  g-th  iteration.  An 
initial  design  defined  by  X  must  be  supplied.  The 
optimization  process  than  proceeds  in  two  steps.  First, 
the  direction  S,  which  improves  the  design,  is  found,  and 
second,  the  scalar  a  ,  is  determined  which  improves  the 
design  as  much  as  possible  when  moving  in  this  direction. 
The  process  is  repeated  until  there  is  no  further  design 
improvement,  indicating  that  this  is  the  optimum  attainable 
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design.  ?or  further  details  see  Ref.  7. 

D.  COPES  AND  SOBEOUTINE  ANAL12 

In  order  to  simplify  the  use  of  CONHIN  and  to  further 
aid  in  the  design  optimization  process  a  Control  Program  For 
Engineering  Synthesis,  COPES,  was  developed  by  Vanderplaats 
[Ref.  7].  COPES  is  the  main  program  (recall  that  COHMIM  is 
written  in  subroutine  form) .  The  user  must  supply  an 
analysis  subroutine  with  the  name  ANALI2,  which  will 
calculate  the  various  parameters.  This  subroutine  has  three 
segments;  INPOT,  EXECUTION,  OUTPUT. 

All  parameters  which  may  be  design  variables,  objective 
functions  or  constraints  are  contained  in  a  single  labeled 
common  block  called  OLOBCa. 

Copes  Terminology 

The  COPES  program  currently  provides  six  specific 
capabilities ; 

1.  Simple  analysis,  just  as  if  COPES  was  not  used. 

2.  Optimization-Minimization  or  maximization  of  one 
calculated  function  with  limits  imposed  on  other 
functions. 

3.  Sensitivity  analysis-  The  effect  of  changing  one  or 
more  design  variables  on  one  or  more  calculated 
functions. 

4.  Two-variable  function  space-Analysis  for  all  specified 
combinations  of  two  design  variables. 

5.  Optimum  sensitivity-  The  same  as  sensitivity  analysis 
except  that,  at  each  step,  the  design  is  optimized  with 
respect  to  the  independent  design  variables. 

6.  Approximate  optimization-  Optimization  using 
approximation  techniques.  Usually  more  efficient  than 
standard  optimization  for  up  to  10  design  variables  or 
if  multiple  optimizations  are  to  be  performed  [Ref.  7]. 
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IV.  DESCaiPTIOa  op  th-e  computer  program 

A.  GENERAL  PROGRAM  ORGANIZATION 

A  functional  bloclc  diagraa  of  the  program  is  presented 
in  Fig.  U.  A  general  description  of  the  subroutines 
contained  in  the  program  is  given  here.  Appendices  A 
through  D  discuss  the  preparation  of  input  data,  list  the 
important  computer  program  nomenclature,  and  list  the 

program. 

B.  SOBROariNES 

1 .  Subroutine  ANALI2 

Subroutine  AMALIZ  organizes  the  basic  anal/sis  used 
in  the  optimization.  It  controls  the  reading  of  the  initial 
design  description  and  calculation  of  the  values  of  the 
objective  function,  constraints,  and  all  other  parameters 
necessary  to  solve  the  problem.  COPES/COMfllN  updates  the 
design  to  minimize/maximize  the  objective  function, 
iterating  until  no  further  improvement  in  the  objective 

function  is  possible  without  violating  one  of  the 
constraints.  COPES/CONMIN  calls  subroutine  ANALIZ  to  obtain 
the  function  value  during  the  optimization. 

2.  Subroutine  INPOT 

This  subroutine  reads  all  input  data  associated  with 
the  bralte  analysis.  Instructions  for  problem  deck 

preparation  are  given  in  appendix  B. 

3.  Subroutine  ISMPR 

This  subroutine  calculates  the  heat  transfer 
constants  such  as  the  thermal  capacity  of  each  node  and  the 
resistance  of  each  node,  determines  the  time  increment  in 
order  to  insure  a  stable  solution,  and  calculates  the  rate 
of  heat  generation.  In  order  to  calculate  the  temperature  of 
each  node,  it  calls  two  subroutines.  From  subroutine  BRAK  it 
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obtains  the  deceleration  needed  to  calculate  the  rate  of 
heat  generated  and  from  subroutine  TEMA  it  obtains  the 
temperature  rise  of  each  node.  Then  it  calculates  the 
temperatures  during  the  time  that  the  brake  is  not  in  use. 
This  subroutine  is  also  capable  of  calculating  the 
temperature  rise  of  a  drum  when  a  constant  rate  of  heat 
dissipation  is  given. 

4.  Subroutine  rsiA 

This  subroutine  calculates  the  temperature  of  each  node. 
As  mentioned  before,  the  heat  is  generated  on  the  inner 
surface,  and  on  the  outer  side  of  the  drum  the  heat  is 
dissipated  by  convection.  The  formulas  used  were  developed 
by  the  finite  difference  method,  and  are  given  in  section 
II-E-2. 

5.  Subroutine  BRAK 

This  subroutine  calculates  the  torque,  actuating  force, 
and  the  friction  moment  of  one  shoe.  It  also  calculates  the 
drum  volume  and  the  deceleration  of  the  machine.  The 
subroutine  takes  into  consideration  a  constant  friction 
coefficient  until  a  temperature  of  90°C  is  reached  and  a 
linear  decrease  in  the  friction  coefficient  for  higher 
temperatures.  More  details  are  given  in  section  II-D-4. 

6.  Subroutine  OUTPUT 

This  subroutine  echos  the  input  data  and  prints  out 
the  thermal  and  mechanical  information  for  the  brake.  An 
example  of  the  output  obtained  from  this  subroutine  is  shown 
in  Table  7  and  Table  2. 
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V.  TEST  PaOBLEM  AMD  HES0LT3 


The  cDaputer  program  was  tested  with  the  data  specified 
in  Table  1.  The  objective  function  which  was  minimized  was 
tna  volume  of  the  drum  material.  Design  variables  were  the 
drum  width,  the  angle  between  the  hinged  pin  and  the  end  of 
the  lining,  the  ratio  of  the  pivot  to  center  of  rotation 
distance  to  drum  radius,  and  the  drum  thickness.  The  side 
constraints  (limits^  on  the  design  variables  were; 


Design 

Lower 

Upper 

yariable 

Bound 

Bound 

1. 

(3)  Width, b 

0.0 

80  mm. 

2. 

(5>  Theta  2, 

1.2  rad. 

2.5  rad. 

3. 

(12) Ratio, Bd 

0.1 

0.  9 

4. 

(19)  Thickness, tk 

40  mm. 

Mo  bound 

The  number  in  parentheses  is  the  location  of  the 

variable  in  the  COasOM  block  in  the  computer  program. 

Constraints  were  imposed  on  the  actuating  force  F,  the 

maximum  temperature,  T  ,  on  the  inner  surface  of  the  drum 

max 

and  stopping  time,  t. 

Constrained  Lower  Upper 

y  ariable  Bound  Bound 

1.  (9|  Force, F  200.0  M-m  2500.0  M-m 

2.  (25|  Time,  t  So  bound  7.00  Sec. 

3.  (!»|  Tempe rature ,T  No  bound  230. 0  “c 

The  vehicle  which  weights  25700.0  Newtons  is  stopped 
four  consecutive  times  from  a  velocity  of  90.0  Km/hr  to 
zero,  with  an  acceleration  period  of  20.0  sec.  between 
stops.  The  values  of  the  design  variables  and  the 
constraints  before  and  after  optimization  are; 
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Before  After 

Jptiaization  Optimization 


Objective 

Function 

Drum  yolurne 

0.754  B-03  m^ 

0.159  E- 

Design 

variables 

Width 

0.03  m 

0.08  m 

Theta  2 

2.10  Rad. 

1.92  Sad 

a/r 

0.75 

0.755 

Thickness 

0.010  m 

0.020  a 

Constraints 

Actuating 

Force 

2815.1  N 

2086.3  N 

Stopping 

time  (last  stop) 

7.04  sec. 

7.00  sec 

Temperature 

after  Last  stop 

348.7  '>C 

229. 200 

Note  that  the  objective  function  increased  as  a  result 
of  optimization.  This  is  because  the  initial  design  violated 
constraints  on  stopping  time  and  maximum  temperature. 

Further  results  are  listed  in  Tables  1  and  2.  In 
addition  to  optimization,  a  sensitivity  analysis  of  the 
design  variables  and  a  two-variabla  function  space  analysis 
for  width  and  thicicness  were  performed.  The  graphical 
results  ace  given  in  Figs.  5  through  17.  The  results  can  be 
summerized  as  follows; 

a.  The  affect  of  changing  the  inside  drum  radius  with  all 
other  design  variables  held  constant; 

As  shown  in  Figs.  5-7,  for  small  inside  drum  radii  the 
drum  temperature  is  very  high.  The  stopping  time  is 
long  and  the  torque  is  low.  Inside  drum  radii  over  130 
mm  give  reasonable  drum  temperature  and  stopping  time, 
for  the  example  considered. 
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b.  As  seen  in  Figs.  8-10,  the  affect  of  changing  the  drao 
width  with  all  other  design  variables  held  constant  is 
the  same  as  described  above. 

c.  The  effect  of  changing  the  drum  thichness  with  all 
other  parameters  held  constant  is;  For  a  drum  thickness 
up  to  6  mm,  the  stopping  tiae  and  drum  temperature  are 
considerably  high.  Over  16  am  thickness,  the  stopping 
tima  re  nains  almost  constant.  For  a  small  thickness  the 
torgue  is  very  low  due  to  the  high  temperatures.  For 
thicknesses  over  20  mm,  the  torgue  remains  about 
constant. 

d.  The  effect  of  changing  the  angle  between  the  hinged  pin 
and  the  end  of  the  lining  is;  For  a  small  02  angle  the 
stopping  tiae  is  very  long  because  the  torque  is  low. 
The  stopping  tima  becomes  raasonable  when  02  >1.8  Rad. 
Obviously  there  is  an  increase  in  the  drum  temperature 
as  02  increases  but  the  ovacall  change  in  temperature 
is  small. 

a.  From  the  two  variable  function  space.  Fig.  17,  it  can 
be  seen  that  the  constant  volume  line  and  the  constant 
temparatura  line  are  almost  parallel,  this  leads  to  the 
conclusion  that  for  the  cycle  taken,  the  drum  is  a  heat 
sink,  and  the  amount  of  heat  dissipated  by  convection 
during  this  cycle  is  small. 
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TEMPESATUaS  AAISS  -  SIMPLIFIED  CALCDLATIpy 

A  siiplified  way  of  finding  the  temperature  rise  of  the 
drum  is  by  using  the  equation; 


Q  =  |cAT 


(40) 


and  setting  Q  equal  to  the  amount  of  heat  generated  using 
Equation  (21)  from  section  II-D-4.  This  equation  is  in 
common  use  in  engineering  texts  (Sea,  for  example.  Refs.  2 
and  8),  The  temperature  rise  calculated  this  way,  is  the 
average  temperature  of  the  drum,  and  not  the  temperature  on 
the  interface,  which  can  be  much  higher  (depending  on  the 
rate  of  heat  generated) .  Extreme  temperature  gradients  cause 
distortion  and  excessive  surface  wear.  Therefore  it  isn't 
always  acceptable  to  use  the  simplified  formula.  Prom 
experience,  it  has  been  found  that  the  surface  wear 
increases  dramatically  as  interface  temperatures  approach 
430  to  533  V  (205  to  260  ’’c)  ,  [Ref.  8]. 

A  comparison  of  the  temperatures  calculated  on  the  inner 
drum  surface,  outer  drum  surface,  and  the  average  drum 
temperature  calculated,  using  equation  (40),  is  given  in 
Fig.  19.  The  graph  shows  the  temperature  rise  for  a  vehicle 
stopped  from  a  velocity  of  90  icm/hr.  Prom  this  graph,  it  can 
be  seen  that  the  drum  temperature,  based  on  equation  (40), 
after  the  vehicle  stopped  is  about  the  average  temperature 
of  the  inner  and  outer  surface  temperatures. 

The  results  show  that  the  drum  will  reach  an  uniform 
temperature  of  about  58  “c,  in  15  sec,  after  the  vehicle  has 
stopped. 

Calculating  the  temperature  with  the  simplified  formula, 
can  lead  to  errors  in  the  time  needed  to  stop  the  vehicle. 
Because  the  temperature  calculated  with  the  simplified 


f^rnula  is  lower  than  the  tenperature  at  the  friction 
interface,  the  calculated  friction  coefficient  is  higher 
than  the  actual  friction  coefficient.  Therefore  the 
calculate!  stopping  tiae  will  be  shorter  than  the  real 
stopping  tiae.  Ail  this  is  true,  provided  the  friction 
material  behaves  as  assumed  in  section  I-D-4. 

Because  high  temperature  is  detrimental  to  both  the 
stopping  ability  and  the  wear  characteristics  of  the  brake, 
it  is  important  that  the  interface  temperature  be  calculated 
with  reasonable  accuracy  in  design.  Fig.  18  clearly  shows 
the  temperature  differences  resulting  from  the  two 
approaches. 

This  difference  in  results  is  compounded  when  the 
simplified  equation  is  used  for  design.  Table  3  presents  the 
design  results  based  on  the  simplified  approach.  This  design 
represents  an  apparent  material  savings  of  27%.  However, 
when  this  optimum  is  analysed  using  the  finite  difference 
heat  transfer  solution,  the  maximum  temperature  is  268.8 ‘’c 
and  the  last  stopping  time  is  7,54  sec.  This  tiae  violates 
the  constraint  by  1.1%.  Perhaps  more  importantly,  the 
temperature  at  the  interface  of  about  269 °C  would  surely 
lead  to  premature  failure.  Therefore  this  design  is  clearly 
too  unconservative  to  be  acceptable. 
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VII.  CONCL'JSION 


In  sammary,  a  nuaerical  optimization  program  is  an 
effective  way  of  finding  a  solution  to  an  engineering 
problem,  provided  reasonable  care  is  used  in  formulating  the 
problem. 


VIII.  aECOMHBMD&TIONS  FOR  FaiDRE  INVESTIGATIONS 

The  study  has  shown  the  feasibility  of  using  numerical 
optimization  in  the  design  of  Internal-Expanding  Rim  Brakes 
with  two  leading  shoes.  Further  studies  on  the  same  design 
may  be  pursued  by  eliminating  some  of  the  restrictions.  For 
example: 

1.  To  add  heat  dissipation  by  radiation. 

2.  To  investigate  drum  temperatures  for  a  drum  with  fins, 

3.  To  take  into  consideration  changes  in  the  surface 
pressure  as  a  function  of  friction  coefficient. 

4.  To  repeat  all  the  calculations  for  a  drum  in  which 
there  is  one  trailing  shoe  and  one  leading  shoe. 

5.  To  add  the  effect  of  centrifugal  forces  for  clutches. 
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LIST  OF  PARAMETEPS 

A  complete  listing  and  description  of  all  variables  used 
in  the  program,  is  not  practical.  The  variables  listed  in 
this  appendix  are  common  to  several  subroutines  of  the 
program  and  will  assist  the  reader  in  a  study  of  the 
program.  The  Global  location  is  the  location  of  the 
parameter  in  the  common  block  called  GLOBCM.  This  common 
block  is  the  means  by  which  Information  is  transfered 
between  the  subroutines  and  the  GDPBS/CONMIM  program. 


Global 

Location 

Fortran 

Name 

Math. 

Symbol 

Definition 

1 

PI 

r 

Inside  drum  radius  (m) 

2 

T 

To 

Torque  of  one  shoe  (N-m) 

3 

NOTH 

b 

Drum  width  (m) 

4 

PRSA 

P 

Pressure  between  lining 
and  drum  (N/m^) 

5 

TETAI 

di 

The  angle  between  the 
hinged  pin  and  the  (Rad.) 
begining  of  the  lining 

6 

XETA2 

02 

The  angle  between  the 
hinged  pin  and  the  end 
of  the  lining  (Rad.) 

7 

PRMNT 

Mf 

Friction  moment  (N-m) 

8 

ANNRT 

an 

Normal  momemt  (N-m) 

9 

ACFPC 

p 

Actuating  force  (N) 

10 

C 

d 

Distance  from  actuating 
force  to  the  hinged  pin(m) 
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11 

Q 

Q 

12 

8D 

a 

13 

CMIO 

14 

HMIO 

15 

AMIO 

u 

16 

SRFC 

17 

BO 

10 

THK 

tk 

19 

OX 

20 

RTIEB 

R 

21 

W 

w 

22 

DCCE 

dc 

23 

TOT 

24 

SCEN 

a 

25 

NWRT 

26 

TIME 

27 

70L 

Vo 

28-32 

TEPL(5) 

T 

33 

HWB 

34 

NWRA 

35 

NWRQ 

36 

NEL 

37 

NS  EG 

38 

?I 

TT 

39 

G 

g 

40 

K 

k 

41 

HCMV 

h 

42 

SPHT 

c 

Heat  generated  (J/sec.) 
Distance  from  pivot  to 
center  of  rotation  (m) 
Cold  friction  coefficient 
Hot  friction  coefficient 
Friction  coefficient  at 
any  temperature 
Drum  surface  area  (m^) 
Outside  drum  radius  (a) 
Drum  thicicness  (m) 
kn  incremental  thickness 
(m) 

Wheel  radius  (»> 

Car's  weight  (N) 
Deceleration  (m/sec/  ) 
Total  time  (sec.) 
Scentricity  (m) 

Write  statement  control 
rime  (sec.) 

Drum  volume  (m^) 
Temperature  at  time  p+1 
(sec.) 

Write  statement  control 
Write  statement  control 
Write  statement  control 
Humber  of  elements 
Humber  of  segments 
Constant 

Sravitational  constant 
Thermal  conductivity 
(J/m-C) 

Convection  heat  coef f . 
(W/m^-C) . 

Specific  heat  (J/kg-^) 
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43 

BHO 

P 

44 

DT&D 

45 

DFTM 

T 

46-50 

TEMP 

T 

51-57 

3ES 

58-63 

TC 

C 

64 

BIO 

3i 

65 

FUR 

?o 

66-72 

NVT 

73-79 

VT 

80 

SELO 

81 

NELT 

82 

TETAA 

0 

a 

83 

ACOF 

84 

TINI 

T 

85 

ZHAN 

t 

86 

NS  HO 

Density  (Kg./nM 
rime  increment  (sec.) 
ilax.  temp,  difference  ( ‘^C) 
Temp,  at  time  p  (sec.) 
Heat  resistance  ( °C/J) 
Heat  capacity  (J/°C) 

Biot  moduli 
Fourier  moduli 
Control  parameter 
Control  parameter 
Number  of  elements+1 
Number  of  elements+2 
The  angle  at  which  the 
pressure  between  the 
lining  and  drum  is 
maximum.  (Bad.) 

Constant 

Initial  temperature ( °C) 
Time  increment  (sec.) 
Number  of  shoes 
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APPEHDIX  B 

IMSTROCTIONS  ?Oa  PROBLEM  DATA  PREPARATION 

Although  the  procedure  is  straight  forward,  preparation 
of  input  data  for  the  program  requires  attenrion.  Errors  are 
easy  to  make  and  difficult  to  locate.  Input  data  is 
dascribed  here  for  the  brake  analysis.  Por  instructions  on 
data  preparation  for  optimization  see  Ref.  7.  Input  data 
should,  in  general,  follow  the  steps  outlined  below.  The  use 
of  the  standard  FORTRAN  Eighty  Columm  Coding  Sheet  is 
recommended.  Integer  constants  must  be  right  justified  in 
the  appropriate  field.  There  are  eight  input  cards,  read  by 
subroutine  INPUT,  to  describe  the  initial  design,  material 
properties  and  constants.  Card  format  is  given  in 
parenthesis  followed  by  specific  instructions  where 
necessary. 

1.  First  Card  <11 0)  -  Duty  cycle  information. 

Cols  1-10  ;  Total  number  of  consecutives  stops 

and  accelerations  (MSEC) 

2.  Second  Card  (I10,3F10.0)  -  Duty  cycle  information. 

a.  Cols  1-10  :  Control  number. 

1  means-deceleration, 

2  means-brake  not  in  use, 

b.  Cols  1 1-20  :  Velocity  at  start  of  deceleration. 

c.  Cols  21-30  :  The  velocity  at  the  end  of  the 

deceleration. 

d.  Cols  31-40  :  The  time  the  brake  is  not  in  use. 

3.  Third  Card  (5110)  -  Thermal  analysis  information. 

a.  Cols  1-10  :  Number  of  nodes  (NEL)  . 

b.  Cols  11-20  ;  An  integer  number  that  controls  the 

amount  of  printout  when  detailed 
output  is  required  during  the 
vehicle  deceleration.  The  amount  of 
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4. 


5. 


c.  Cols  21-30 


a.  cols  31-40 


e.  cols  41-50 

Fourth  Card  (7 

a.  Cols  1-10 

b.  cols  11-20 

c.  Cols  21-30 

d.  Cols  31-40 


9.  Cols  41-50 

f.  Cols  51-60 

g.  Cols  61-70 
Fifth  Card  (7F 
a.  Cols  1-10 


lines  written,  depends  on  the 
stopping  time  and  time  increment. 
(NHR)  . 

:  An  integer  number  that  controls 
the  amount  of  printout  when 
detailed  output  of  the  temperatures 
is  required  during  the  period  that 
the  vehicle  is  not  in  use.  The 
amount  of  lines  written  depends  on 
the  period  length  that  the  brakes 
are  not  used  (NWRA)  . 

:  An  integer  number  that  controls  the 
amount  of  printout  when  detailed 
output  of  the  temperatures  are 
required  during  the  period  of 
constant  heat  generation.  (NMRQ) . 

;  Humber  of  braking  shoes  in  the 
machine. 

F10.0)  -  Brake  dimensions. 

:  Inside  drum  radius  (HI) . 

:  Drum  width  (WDTH)  . 

:  Drum  thickness  (THK) . 

:  Ratio  of  distance  from  pivot  to 
center  of  rotation  and  inside 
radius  (RD)  . 

:  Drum  density  (RHO) . 

:  Angle  between  binged  pin  and  the 
begining  of  the  lining  (TSTA1)  . 

:  Angle  between  hinged  pin  and  the 
end  of  the  lining  (TETA2)  . 

10.0)  -  Thermal  and  friction 
information. 

;  Heat  conduction  coefficient  (K) . 
(real  number) . 
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b.  Cols  11-20  :  Heat  convection  coefficient  (HCNV) . 

c.  Cols  21-30  :  Specific  heat  of  the  drum  (SPHT) . 

d.  Cols  31-40  :  Hax.  temperature  difference  between 

cold  friction  coefficient  and  hot 
friction  coefficient  (OFTM) . 

e.  Cols  41-50  :  Cold  friction  coefficient  (CMI7) . 

f.  Cols  51-60  :  Hot  friction  coefficient  (HMIV)  . 

g.  Cols  61-70  :  Initial  temperature  (TINI) . 

6.  Sixth  Card  (2P10.0)  -  Machine  information. 

a.  Cols  1-10  :  Vehicles  weight  (H)  . 

b.  Cols  11-20  :  Wheel  radius  (RTIES)  . 

7.  Seventh  Card  (5P10.0)  -  Analysis  constants. 

a.  Cols  1-10  :  Maximum  pressure  between  lining 

and  drum  (P5SA)  . 

b.  Cols  11-20  :  Constant  3.1415927 

c.  Cols  21-30  :  Gravitational  constant  (G) . 

d.  Cols  31-40  ;  Increment  of  time  (ZHAM). 

e.  Cols  41-50  :  The  angle  of  maxiouai  pressure 

(TETAA)  . 

8.  Eight  Card  (110.0)  -  Print  control. 

Cols  1-10  :  An  integer  number  can  be  zero  or  1. 

If  zero  (or  a  blank  card)  -  only 
the  final  results  are  printed. 

If  1-  the  temperature  at  time 
increments  are  printed. 
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APPENDIX  E 


0)1 


Clutch  or  Brake 


U)  2 


Fig.  1  Dynamic  Representation  of  a  Brake  or  Clutch 
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Inner  Surface 


Fig.  3  Finite  Difference  Model 
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Fig.  4  Block  Diagram  of  the  Program 
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Fig.  5  Maximum  Inside  Drum  Temp.  Vs.  Inside  Drum  Radius 


Fig.  6  Stopping  Time  Vs.  Inside  Drum  Radius 


Fig.  7  Torque  Vs.  Inside  Drum  Radius 


Vs.  Drum  Width 


Fig.  9  Stopping  Time  Vs.  Drum  Width 


Torque  Vs.  Drum  Width 
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Fig.  11  Maximum  Inside  Drum  Temp.  Vs.  Drum  Thickness 
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Fig.  12  Stopping  Time  Vs.  Drum  Thickness 


Fig.  14  Maximum  Inside  Drum  Temp.  Vs.  Theta 
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